Continuity in two dimensions for a very degenerate elliptic equation 
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' Abstract: An elliptic equation V • (-F(Vn)) = / whose ellipticity strongly degenerates for small 

. values of Vu (say, F = on -8(0, 1)) is considered. The aim is to prove regularity for F{Vu). The 
paper proves a continuity result in dimension two and presents some applications. 

<^ . 1 Introduction 

■ 

+^ . We consider in this paper an elliptic equation of the form 

V-{F{Vu)) = f, (1.1) 

where F is the gradient of a convex function H : ^ R. This is one of the most classical 
^ , non- linear elliptic equations, which arises as an optimality condition for the minimization of the 

QQ ■ functional / H{Vu) + J fu. When H{z) = |zp or H{z) = \z\p we get the usual Laplace and 
. p— Laplace equations, respectively, which have been investigated intensively in the literature and 

(N 



have provided a lot of regularity results on the solution u (which is a priori supposed to belong to 
or W^'P only) according to the regularity of /. Most of the results have been extended to the 
■ case of variable coefficients (see the book by Gilbarg and Trudinger [6] for a compendium of the 
0^ ! theory of elliptic regularity and the original paper by De Giorgi, [3], where he proved a key Holder 
regularity result) or of different functions H, which share anyway some properties of the square or 
of the p—th power. This latter case is much more difficult than that of the square (which gives a 
^ , linear equation), mainly because of the degeneracy of D^H near z = 0. Actually, lower bounds on 
I D^H are often useful to give estimates on the norms of the solution in terms of the norms of / 
and if D^H is allowed to tend to zero for small values for Vn (which is the case fro p > 2), extra 
difficulties arise. Yet, this is - roughly speaking - compensated by the fact that the smallness of the 
gradient already provides some regularity estimates and some results are still possible. 

This is what is usually done in elliptic regularity but the equation we want to look here is even 
worse than the Laplace equation. Our main example is given by the function 

H^g){z) = ^i\z\-l)l, q>l (1.2) 
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which vanishes, together with its Hessian, on the whole ball Bi. This means that all the values of 
Vu which are smaller than 1 are a source of problems. Obviously, the regularity result that one may 
expect to prove must concern F(Vu) instead of u itself: it is evident that the homogeneous equation 
(the same, with / = 0) would be solved by any 1— Lipschitz function and that nothing more may 
be said on u. On the contrary, F(Vu) is reasonably more regular, since either it vanishes or we are 
in a zone where |Vu| > 1 and the equation is more elliptic. But gluing the two zones, which are not 
open sets with nice boundaries, is not trivial. 

Up to our knowledge, the first need for a regularity result on F(Vu) in the case of a degenerate 
function F of this kind has been encountered in [Ij for a variational problem linked to traffic 
congestion. In this problem, i.e. 

min jy + ^\cr\'^ : a £ L'^{n), V ■ a = f, cj • n = o| (1.3) 

the optimal a equals -F(Vti) where n is a solution of (jl.ip with F = V-ff(g), Neumann boundary 
conditions and ^ + | = 1. For the sake of the applications of [T] boundedness and Sobolev regularity 
were all that was needed on a, and they are proven for functions H of the form given in (jl.2p . 
In particular it is proven that a = F(Vu) belongs to n in any dimension d. Continuity is 
not addressed, even if it would be interesting to obtain. We will explain in Section 6 the possible 
applications of a continuity result, both in the framework of traffic congestion and of the variational 
problem (jl.Sp in general. 

Here we consider the more general case of a function H whose Hessian is bounded from below 
outside any ball Bi^s and we prove a continuity result which is valid in two dimensions, provided 
that we already know that F{\/u) belongs to n L°°. This is exactly the case for H = H^^^y Yet, 
the assumptions on / are more standard in this framework (/ € i^"*"^), while the Sobolev results of 
[1] asked for Sobolev regularity of / itself, which is not that natural in elliptic regularity. 

The proof of the present paper is two-dimensional because it follows a method developed by 
DiBenedetto and Vespri for different equations, which is based on the following idea: using the 
equation, one can prove that if the oscillation is not significantly reduced when passing from Bn to 
-Beq/j, then the contribution of the crown B^ \ -BeoR to the Dirichlet energy (which is finite because 
of the assumption) is at least a certain value. This value, in the case of d = 2, scales in a way 
so that it does not depend on R. This is used so as to prove a decay for the oscillation and get in 
the end a logarithmic modulus of continuity. 

In our case we first prove results for functions of the type {du/dxi — (1 + (5))+ and then deduce 
some results at the limit as 5 ^ 0. This gives continuity of all the expression ^(Vu) for all the 
functions g vanishing on Bi. Typically, this includes -F(Vu) itself. 

Notice that all the results we present are local and that we did not try to improve them up to 
the boundary. 

The strategy followed by the paper and its plan will be detailed in Section 2 below. 
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2 The equation, the assumptions and the strategy 



Let : M°' ^ M be a convex function, satisfying the key assumption that 

for all 5 > there exists cs > : D'^H{z) > csid for all z such that \z\ > 1 + 6, (2-1) 

where Id is the identity matrix. The function H may be for simplicity supposed C^, (or Cj^l). 
For notational purposes, as we already did in the introduction, we will denote by F, so that 
: R'^ — > M*^ is a vector function satisfying some monotonicity assumptions. We will deal with 
solutions of the following equation: 

V • {F{Vu)) = /, 

being / a given datum in a suitable space. We will see that for most of the results we prove L-^+^ 
will be sufficient, while for some of those that we recall we would need / in a Sobolev space. 
Notice that this is the equation which is satisfied by the minimizers of 



mm 



H{Vu) + I fi 



(we do not precise boundary conditions since we are only interested in local results). The function H 
is allowed, for instance, to vanish on the whole Bi, as it is the case for the typical example H = H(^gy 
Hence, H is in general not strictly convex and no uniqueness is valid for the minimization problem 
above. Yet, it is easy to see that F{Vu) is the same for all the minimizers. This vector field will be 
the main object of investigation. By the way, F^Vu) is also the solution of the dual problem 



min J H*{a) : V ■ a = f. 



Strategy. For a solution u of our problem, we write down the equation which is satisfied by a 
partial derivative v = du/dxi. 

To do so, we are tacitly assuming that u is regular (at least -ff^), which is not at all obvious with 
the kind of equation we have. Hence the idea is: approximate the equation by truly elliptic equations 
(for instance changing the function H into a function H^, still satisfying (|2.ip . but satisfying also 
C~ < D^H^ < C^, for some non- uniform constant C^, and fixing suitable boundary values: see for 
instance [l]); then write estimates on the solutions which will only depend on cs and not on the 
constants C^. This estimate will pass to the limit as e — > 0. The functions will converge up to 
subsequences to a solution of the limit equation, which has no uniqueness, but the quantity F{Vu) 
is the same for all the solutions. Hence, estimates on F^Vu) will be true for the non-approximated 
equation as well. 

In the following, we will not precise any more the approximation, but ,with this approximated 
spirit in mind, let us differentiate the equation and get 

V • (a(x) • Vv) = /', 
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being a{x) = D'^H(Vu{x)) the Hessian matrix of H and /' the partial derivative of /. This is a 
Unear ehiptic equation with the following property: a is not uniformly elliptic and can also vanish, 
but on the points where \v{x)\ > 1 + 5 one has a fortiori |Vtt(x)| > 1 + 5 and hence a{x) > csld- 

The main object of our preliminary analysis will be the function vs = hi^s{v) (we define /ifc(t) := 
{t — k)+). Since is a Lipschitz and convex function, vs is a subsolution of 

V • (a(x) • Vvs) = f'lv>i+5- 

This new equation may be assumed to be uniformly elliptic, since the values of the matrix a 
on the points where vs = are not important. Yet, vs is only a subsolution and this will not be 
sufficient to establish all the results we need for regularity. At a given point, we will also need to 
use the equation which is satisfied by v itself (which is not uniformly elliptic). 

In Section 3 we will prove a continuity result for vs, in dimension 2 and under some a priori 
assumption on vs itself. We will need to suppose that vs is a function and that it is bounded. 
Notice that, due to the fact that we tacitly regularize, when we say that vs or / belong to a certain 
functional space, we actually mean that the estimates we prove will only depend on the norm in this 
space, and hence will be inherited by the solutions corresponding to non-regular data, provided we 
stay in the same space. Also the continuity result is given in the same spirit, in the sense that we 
will prove a quantified continuity (in our case, the modulus of continuity will be uj{R) = C\ lni?|~^/^ 
and the constants will depend on q). The techniques follow a strategy by DiBenedetto and Vespri 
(see [5]) which amounts at proving some oscillation decay from a ball to a smaller one under some 
conditions. These conditions are such that violating them implies that the solution accumulates a 
fixed amount of Dirichlet Energy in the crown between the two balls and the assumption allows 
to conclude. 

Section 3 aims at being as general as possible and some of the lemmas will be stated in M*^ instead 

of 

In Section 4 we will see how to link vs to -F(Vm). In particular this will allow to translate Sobolev 
results on F{Vu) into corresponding results for vs, and also to state that vs as well does not depend 
on the non-unique solution u. 

In Section 5 we will conclude wider continuity results on functions of Vu, as a consequence of 
the results on vs- 

Section 6 will present some interesting applications, in cases when we do actually know that 
i^(Vit) is Sobolev and bounded. This is for instance the case for H = (see [1]). The 

regularity mainly depends on structure assumptions on the function H and may fail to be valid in 
general, while the one seems to be easier, since it only depends on the fact that functions of the 
kind {du/dxi — (1 + 5))+ are subsolution of a uniformly elliptic equation. 

3 Continuity for vs in dimension 2 

The section will be opened by some very classical lemmas which have proven to be useful several 
times in elliptic regularity. For brevity, we omit the proofs and we refer to the book by DiBenedetto 
[1]. We also refer to the original work by De Giorgi [3]. 
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Lemma 1. // (l^)n is a sequence of non-negative numbers satisfying 
then Yji ^ 0. 

Lemma 2. The following Poincare-type inequality holds on any domain J7 in with a universal 
constant (depending on the dimension d only): if u E Hq{0,), then 



[ < C|{u/0}|2/'^ [ |Vn|2. 
Jq Jn 



If the previous inequality comes from an idea by De Giorgi, the following is very very standard: 

Lemma 3. // a shape Qi is fixed, then for any rescaled open set fi/j = RQi C M*^ the following 
Poincare-type inequality holds with a constant C depending on and on the ratio X: ifu G H^{^ji) 
and \{u = 0}| > Airi^l, then 

Lemma 4. Suppose f G Lp''^^ , and let v he a suhsolution of 

V • {a{x) ■ Vv) = f'l^>o, 

where the matrix a satisfies cid < a < CId- Then there exists a constant vq such that, ifO<v<M 
in Br and 

\{v > {l-a)M}nBR\ < uoIBrI 

£ 

then V < (1 — q/2)M in Br/2j provided M > R^+^ . The constant vq only depends on the dimension, 
on a, on e and on the bounds and on the ellipticity of a (i.e. on c and C). 

Proof Set 

A;„ = M(^1---— j, R^ = - + —, Ar, = {v>kn}f^BR^, 

let rjn be a cutoff function which vanishes outside Br^^ and equals 1 on Br^_^_^ , with gradient bounds 
|V?7„| < CR~^2'^. The function 0„ := (t; — is a subsolution of 

V- (a(x)- V<^n) = /%>ifc„, 

and we can test the equation against ^n'Vn- We get 

c j \y(t>n?rii < J < aV0„, V</)„ >vl<C J |V0„||V?7n I + J f'lv>K ■ 

The first term in the right hand side may be estimated by | / |V(/)nP?7^ + ^ / ^nl^^nP; in the 
second one can get rid of Iv>k„ because anyway (pn vanishes on {v < and integrate by parts, so 
as to obtain 
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Putting together the estimates one gets 



j \V^n\'r]l <c{^j <Pl\Vr,r,\'' + I frii^ . 
We then apply the estimate of Lemma [2] to the function u = (pniln and get 



In the last expression, we may estimate (j) with M (actually M(l — a) would be sufficient), |Vr/„| 
with 02"^ /R, and apply Holder inequality to / and Ia„- On the other hand, the values of 4>'^rj^ may 
be estimated from below by {kn+i — kn)lA„+i- As a consequence, we get 



^n+l ~ ^n)^|^n+l| ^ C 



p2+E 



2 

2+e 



We define y„ as (i-e., up to possible bounded factors, the ratio between the area of An and 

the ball itself). We have 



Replacing Yn with Yn'^^ (since Yn < 4), and supposing M > R'^+' ^ one gets an estimate of the form 



y„+i < C{d,C,c,\\f\\L2+.,a)Y, 



16" 



If one supposes Yi < i'q, where vq is given by the assumptions of Lemma [H we get y„ — > 0. This 
means that the measure \{v > {1 — a/2)M} n -B_r/2| is zero, which is the thesis. □ 

Lemma 5. Let < u < M be a function in H^{Bji), assume that \{u > 3/4M} n Bji\ > uqIBjiI, 
and set eq = \J vq/2, then one of the following alternatives happens: 

• either the Dirichlet energy of u in Br \ Bi^^^ji is large: /g^^^ ^ |Vu|^ > cM^, 

• or there exists a radius s G [sqR, R] such that u > |M on dBg. 

Proof. Thanks to the value we have chosen for the set {u > |Af} n {Br \ B^^^r) must include a 
large part of Br \ B^^r and in particular we have|{ti > 3/4M} n {Br \ i^^Q/j)! > ^\Br\. 

Let us define a subset X of [eqR^R] through X = {s ^ [eqR.R] : 3x G dBg : u{x) > |M}. 
We obviously have 

2^R\X\ > 27rsds > \{u > ^M} n {Br \ B,,r)\ > '^\Br\ 
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which gives \X\ > ^R. Now, either the second alternative is verified or for every s G X there exists 
a point X on the circle dBg such that u{x) < |-/Vf . On each one of these circles we would get an 
oscillation of at least M/8 and hence 



M 

— < 



/ \Vu\dn^ < (27rs)i/2 f f iVupd^i 

JdBs \JdBs 



1/2 



which implies 



and hence 



|V^z|W> 



B{xo,R)\B{xo,£oR) 



|Vup > ds \Vu\'^dn^ > \X\ > CuoM^. □ 

Jx JdB, USttR 



From now on, we will denote by eo the constant defined in Lemma [5] obtained for q = 1/4 and 
c = cs (i.e. eo = ^Jv^jl where is obtained from Lemma H]). 

Lemma 6. Let u be in H^{Bji\Bi;gfj) and assume that both the measures \{u > aM}n{Bji\BsQji)\ 
and \{u < bAI} n \ i^eg/j)! are larger than rjR? , then 

/ |Vm|^ > c(r/,eo,a,6)M2. 

Jbb\b,^r 

Proof. Consider the function (j) = {u — bM)^ f\ {a — b)M . It vanishes on a set whose measure is at 
least rjR? and is maximal and equal to (a — b)M on a set whose measure is, again, at least rjR?. 
This implies, thanks to Lemma [3] applied to 17 = i?i \ iJ^o, the inequality 

7]R^{b-afM^ < [ (P"^ <Ci7],eo)R^ [ \V(p\'^ < C{rj,£o)R^ [ |Vu|2. 

JBr\B,^r JBr\B,^^r JBr\B,^r 

The last inequality comes from the fact that </) is a 1— Lipschitz function of u and the statement is 
obtained. □ 

Lemma 7. Let Bn and B^^ji be two concentric balls in 0, and vs = hi^s(uxi) as defined in Section 
2. Then one of the three followings alternatives happens 

• either osc{vs,B so r) < | osc(t>5, 

• or the energy of vs in Bji \ Bfr^^n is large: 

\Vvs\'^ >c[osc{vs,BR)f, 

Br\B,^r 

or we have osc{vs, Bji) < R^ . 
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Proof. Up to a translation, suppose that < < M in B^, where M = osc{vs, Bji). In the ball 
Br the function vs (actually it is vs — inf^^ w^) is a subsolution of 

V • (a(x) • Vv) = /7,>o, 

where the matrix a satisfies a > csl at least whenever v > 0. Hence, one can replace a with csl on 
the set {vs = 0} and suppose that the equation is uniformly elliptic. We start by supposing that 
the inequality M > i?2+e ig satisfied, otherwise the second alternative is realized. Assuming this 
estimate, one can apply LemmaSlwith a = 1/4. We get that either \{vs > 3/4M} n Br] > vo\Bji\ 
or the oscillation on Br/2 is reduced by a factor of 7/8 and the first alternative is realized. 

Suppose now that \{vs > 3/4M} n Br\ > vqIBrI. Lemma [5] proves that in this case either 
the third alternative is realized or there exists at least a radius such that vg stays high on a circle 
between Br and B^^r. 

Hence, suppose now that at least a value s G [sqRjR] is such that vs > |M on dBs- In such a 
case, consider the function w given by 



(vs - |M)_ a iixG Bs, 

iix(^Bs. 



To get information on the behavior of w one needs to look back at the equation satisfied by the 
partial derivative Ux, instead of the one which has vs = [ux — (l + (5))+asa subsolution. We know 
that Ux satisfies 

V • (a(x) • Vux) = /', 

which is a non-uniformly elliptic equation, and we test it against the test function Wn = {w — kn)+, 
where the numbers kn are given by = (1 — 2"^"+^))-^, for n > 4. We get 

f d 

< aVUx,VWn >= - / /t^ Wn 



dxi 

and we notice that, where Vwn 7^ 0, one has Vwn = — Vu^, and also > j^M > 0, which implies 
ffl ^ Csl- Hence we have 



cs 



J \VWn\'^ < J <aVWn,VWn >< J f\VWn\ <J J + ^ / /^^V«;„^0- 

This implies, by using also Lemma [2] 

jwl< C\An\ j \VWn\^ < C|^„|||/||i2+.|^n|^, 

where An = {w > kn}- With the same kind of estimates as in Lemma IU one gets 

\kn+l - Vn+ll < C\An\^+^\\f\\l2+. 
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and, setting again Yn = AnR ^ and using /cn+i — kn = M2 this gives 

The last inequaHty comes from the assumption M > R^+' . Again, if one supposes Yq <vi (where vi 
is given by Lemma [U as in Lemma HI , one gets 1^ ^ 0, which means that w < ^ on Bs, provided 
\{w > %}\ < u^\Bn\. 

This means that there are two remainding cases: either the condition on the measure is fullfilled, 
or not. In the first case the oscillation has gone down to |M on Bg and hence we can say that the 
oscillation on -Beq/j is less than 7/8 of that on Br and the first alternative is realized. 

In the second case, we may estimate again the energy thanks to Lemma [6l In this last case it is 
the third alternative which is realized. □ 

The conclusion of all these lemmas is a continuity result on vg- 

Theorem 8. Let u be a solution of V • {F(Vu)) = f , with f £ L^"'"'' in dimension d = 2, and 
suppose that vs = hi^siuxi)+ S H^nL°° and that D'^H{Wu) < CId (which is always true if H G 
and Vn G L"^). Then vs is continuous on Q, and its modulus of continuity gives \vs{x) — vs{y)\ < 
C (In |x — for any x,y C ^l, the constant C depending on ||/||l2+£:, ||f(5||L°°) ll^^lli^i; 

cs, C, d{VlQ,d^). 

Proof. Take a sequence of balls Bn included in fi, of the form Bn = i?(xo, eo-Ro)i (-^o < diyt^^d^) 
for n = 1, . . . , A^. We call R the last radius, i.e. R = RqEq , and hence 

^ ^ ln(i?/i?o) 



Ineo 

Each time we pass from i?„ to i?n+i either the oscillation is diminished by a factor 7/8, or the 
oscillation was already smaller than a power (3 = e/(2 + e) of the radius, or the energy in the crown 
is controlled from below. Let us denote by j, k and h the number of times the three situations 
happen. We have j -\- k + h = N . If we call M as before the oscillation of v^ on B^f we have the 
three following informations: 

• M{^y < ||f5||L°=, since the oscillation on Bq cannot be more than the maximal value of 

vs, which is positive, and it has diminished j times; we can also write My < ||v<5||l°°) since 
(8)i = (1 + ly > 1 + 2 > 2. 

• M < (RqEq)^ < (RoEq)^, where i is the last index of the second type (and hence one can say 
that Sq < Eq since i > k); here we can write M < C{Ro, P,eo)/k, thanks to the elementary 
inequality < C{£)/t, which is valid for e <1 and t > 0. 

• cM'^H < ll'yall^i, since on each crown Bi \ -Bj+i where the third alternative is realized we had 
a contribution proportional to the squared oscillation (and hence at least M^) to the Dirichlet 
energy of vs, 
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This implies the following bounds on M: 

/ • I 7| 1^-5111,°° \\v5\\]ii C{Ro,I3,£q)\ Tiiax.{l\\vs\\L^;\\vs\\]ji; C{RQ,(3,eo)} 
M < mm < — — : ; — ^ > < 



At least one of the three indices j, k and h is larger than N/S, and hence one gets 



M < C{\\vs\\L^,\H\\m,Ro,s:o,P)^=^ = C(||^^<5||l->, ||^^5||^/i,-Ro,/?,eo)^' 



N \/\lnR\ 

which gives a logarithmic modulus of continuity for vs- This constant degenerates as the ellipticity 

constant cg goes to zero, since Eq = eo{cs,C, f3). The result is local and becomes worse when 
approaching the boundary, due to the dependence on Rq (which can be taken as large as d(i^o, dO,)). 

□ 



4 From F(Vn) to vs 

In this section we show a useful and clarifying lemma, which aims at giving the relation between 
F(Vu) and vs- 

Lemma 9. For any a G M*^ consider all the points z such that F{z) = a and set 

75(a) := (zi - (1 + 5))+ = hi+5{zi), 

where zi denotes the first component of z. 

The function 75 is well-defined (i.e. if several z satisfy F{z) = a the value of hl-^-s{zl) is the 
same for all of them), and it is Lipschitz continuous with a Lipschitz constant that does not exceed 

Proof. Consider points a,b,z,w G so that F{z) = a and F{w) = b. Suppose zi,wi > 1 + S: in 
this case one has, from F = V-ff and from the quantified convex behavior of H on the half-space 
{^1 > 1 + S}, 

< F{z) — F{w), z — w >> cs\z — w\'^, 

which implies |z — w| < cj^\a — b\ (this implies in particular uniqueness of z and w). Then one can 
compose with hi+s, which is 1— Lipschitz, and get 

\hi+s{zi) - hi+s{wi)\ < \zi - wi\ <\z-w\< cj^\a - b\. 

In the case zi,wi < 1 + 5 one trivially has hi^sizi) = hi^s{wi) = 0. The only remaining case is 
(by simmetry) zi < 1 + 5 and wi > 1 + 5. In this case there is a unique point z on the line through 
z and w lying on the hyper plane {zi = 1 + 5}. Hence one may write 

< F{z) — F{w), z — w > > cs\z — w\'^, 
< F{z) - F{z),z- z> > 0, 
<F{z) - F{z),z-w> > 0. 
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The second inequality comes form simple convexity of H and the last one from the second, since 
z — z and z — w share the same direction. Hence, summing up the first and the third, one gets 

< F{z) — F{w), z — w >> cs\z — 

which gives 

wi — (1 + 5) < \z — w\ < cj^\a — b\, 
and, again, the Lipschitz result. □ 

As we briefly addressed in Section 2, this lemma implies that vs is the composition of F(Vti) with 
a Lipschitz function (the Lipschitz constant only depending on S). This is useful in several contexts. 
First of all, we need vs being a function (not necessarily Lipschitz) of F(Vn) if we tacitly approximate 
the equation, since we already underlined that the solutions of the approximated problem may 
converge to a specific solution u of the limit one, but F(Vu) is the same for all solutions. Thanks 
to this lemma, the same is true for vg- Then, this lemma also proves that has the same Sobolev 
regularity of F(Vn), and that its norm only depends on that of F{Vu) and on cs- In particular, 
it does not depend on the direction that we choose for differentiating u (that we conventionally 
called xi but that could have been any other direction). This will be useful in the next section, 
where we will deduce results on |V?x| starting from vg (using the arbitrariness of the direction and 
the equicontinuity of all the functions vs)- 



5 Continuity of functions of Vu 

We have proven quantified continuity results for vs = hiJ^s{du / dxi) . Just changing the variable one 
differentiate with respect to, one gets the same results for every function of the form ^1+5 (Vu ■ e), 
for |e| = 1. Since these functions all share the same modulus of continuity, taking the supremum 
over e one has the same result for /ii+5(|Vu|) = (|Vu| — (1 + 5))+- 

Finally, if one lets (5^0, one gets continuity for (|Vtt| — 1)+ and for (Vti-e— 1)+ as uniform limits 
of continuous functions. The modulus of continuity depends on the dependence of the modulus of 
continuity previously obtained with respect to 5. We did not explicitly accounted for these continuity 
moduli but it could have been possible. Anyway, one can say that the continuity of (|Vm| — 1)+ and 
for (V« • e — 1)+ only depends on the constants {cs)s>o-, on the upper bounds on the Hessian of H 
(in the applications these usually depend on the maximum of |Vu|) and on the L^+^ norm of /. All 
the results are local. 

By composition, it is not difficult to prove continuity for functions of the form ^'dViij), for any 
real continuous function g satisfying g{t) = for t < 1. 

To allow for more general functions and in particular for vector functions of Vu, we present the 
following lemma: 

Lemma 10. Let z : ^0 ^ he a hounded function such that the family of functions h\{z ■ e) is 

equicontinuous for e € M*^ with \e\ = 1. Let g : -^W^ be a continuous function with g = on Bi. 
Then g{z) is continuous on Qq and its modulus of continuity only depends on that of g, on \ \z\\l°° 
and on the common modulus of continuity of the functions hi {z ■ e). 
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Proof. First remark that, as we already pointed out, then the function (jzl — 1)_|_ as well is continuous, 
with the same modulus of continuity, as a supremum of a family of equicontinuous functions. 

Now, fix e > and let 6 = 6{e) > be the corresponding 5 given by the equicontinuity of all 
the functions hi{z ■ e) and ^i(|z|). Take two points x,y £ ^Iq with — y| < 6 and set a = z{x) and 
b = z{y). 

We distinguish two cases: either |a| < 1 + 2e or not. If yes, then we have g{a) < u;(2e), where 
uj : M"*" is the modulus of continuity of g, and 

|(|a| - 1)+ - - 1)+| < 8, with (|a| - 1)+ < 2e, 

which implies (|6|-1)+ < 3e and |6| < l + 3e and g{b) < uj{3e). Finally, \g{a)-g{b)\ < uj{2e)+uj{3e). 

If not, then we have (|a| — 1)+ = \a\ — l > 2e and | (|a| — 1)+ — — 1)+| < e implies — 1)+ > e 
and (|6| — 1)+ = |6| — 1. This gives ||a| — < e. In this case we want to estimate |a — 6|: we use 
the continuity of (z • e — 1)+ with e = a/\a\. We get 



1 



< e 



since, again, we have |a| — 1 > 2e, then (6 ■ 



\a\ 



\a\ 



1 > e and hence 



< e, and thus 



6 • —7 > lal 

a 



Then we write 

|2 I |2 



a 



b\^ = \a\^ + \b\^ - 2a ■ b < ja^ + |6r - 2|a|(|a| - e) = 2|a|e + |&|^ 



2\a\e + {\b\ + \a\ 



\a\) < e(2|a| + \a\ + \b\) < e {4\a\ + e) 



Since z was supposed to be bounded, we can assume \a\ < M and hence \a — 6p < Ce. In this case 
we get \g{a) - g{b)\ < uj{VCe). 

In all the cases we get continuity of g o z. □ 

The following theorem is just a corollary of Lemma [TOl Lemma [9] and Theorem [51 

Theorem 11. Let u be a solution of V ■ (-F(Vn)) = /, in dimension d = 2, where F = VH and 
H is a convex function satisfying the assumption of Section 2, with D^H bounded on bounded sets. 
Suppose f G L'^^^ and that F{Vu) £ n L°°. Then g{Vu) is continuous on Q for any continuous 
function g : — > M such that g = on Bi. 



6 Applications 

As a particular consequence of the previous section, in the case H = (for q > 2), one has 
continuity for F^Vu) itself, since the function F vanishes on the ball Bi and all the assumptions for 
applying the previous theorems are verified (in particular, pj proves and L°° regularity). 
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As we said, a = F{SJu) optimizes J H*{(t) under a divergence constraint, and in this case we 
have HI^){(t) = kl + l\a\P. We summarize our results - under the variational formulation - in the 
following theorem. 

Theorem 12. Fix p < 2. The unique minimizer a of 

min \a\ + ^\a\P : a € LP{n), V ■ a = f, cj • n = o| (6.1) 

is a continuous function in the interior of Q, provided f G Ty^'^(J7). 

Proof. All the ingredients for this result have already been proved. We know that the optimal 
solution a equals F(Vu) where F = Vi?(q) (with ^ + p = 1 ^^'^ Q ^ 2)and u is the solution of 
V • (F(Vu)) = /. Thanks to what proven in [T], the vector field F^Vu) is both bounded and H^. 
This last regularity result requires, for the proof which is performed in [1], the Sobolev assumption 
on /. The result is presented in [1] for a Lipschitz / but it should not be difficult to adapt to 
/ E On the contrary it does not seem easy to relax the assumption for Sobolev regularity. 

Then we use Lemma [9] to say that vs S H^{Q) as well, and that the norm of vs does not 
exceed c_^^||F(V)||/^i (in particular, for fixed 6, it is uniform with respect to the direction e). 

This allows to use the results of Sections 3 and 4 for getting first the continuity of vg and then 
of FiX/u), since in this case F vanishes on the ball Bi. 

^{q) ( V"") comes from H(^g^ G C/j^ 



Boundedness of D'^H/„\ (Vu) comes from H/„\ G and boundedness of Vu. □ 



This result is interesting in itself as a regularity result for a variational problem under divergence 
constraints. 

We also briefly present the consequences that our continuity proof has in the domain of continuous 
traffic congestion, as in [2] and later in [T]. 

The latter paper mainly aims at exploring the equivalence between Problem (j6.ip and the min- 
imization of a total traffic intensity J H* {iQ(x))dx where iqix) denotes the quantity of traffic at x 
associated to a probability distribution Q on the space of possible paths (the continuous framework 
is described through the choice of all Lipschitz curves in 17 as a set of admissible paths, instead of 
taking the paths on a given network). For modelling reasons, the choice of H = H(^q^ or in general 
of a function H such that (-ff*)'(0) > makes more sense than a simple power, since (i?*)'(0) 
represents the cost for passing through a point with no traffic, which should be low but non-zero. 

It turns out from the analysis in [Ij that one can build a measure Q which is optimal for the 
traffic intensity problem starting from an optimal vector field a for (16.1jl . We then follow the integral 
curves of a vector field 

^^^'"^^ (l-t)/+(x)+t/-(x)' 

where and /~ are the positive and negative parts of /, respectively. This would work fine 
provided the vector field is Lipschitz continuous, so that its integral curves are well-defined as the 
solutions of a standard Cauchy problem. Since a may not always be assumed to be Lipschitz, due to 
the degeneracy of the equation that one gets with H = ifj-g) , then the efforts in [T] were devoted to 
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the proof of Sobolev regularity so as to apply DiPerna-Lions theory for weakly differentiable vector 
fields. 

Yet, in the case where a turns out to be continuous (and the present paper proves such a result in 
dimension two), its integral curves may be at least interpreted as classical solution of an ODE with 
non-uniqueness and their regularity could be improved (solutions of a first-order Cauchy problem 
with continuous data are C^). 

Besides this first application, there is another interesting consequence in traffic congestion which 
is more linked to the topics of [2]: it is proven in such a paper that optimal traffic distributions Q 
satisfy a Wardrop equilibrium principium (which says that every path which is actually followed by 
somebody must be actually geodesies with respect to a metric taking into account traffic congestion: 
it is an equilibrium since the metric itself depends on Q and Q must be concentrated on zq— optimal 
curves). Yet, the metric one obtains is a function g[iQ{x)) which is a priori neither continuous nor 
semi-continuous, since iq is only U*. The efforts of a large part of ^ aim at giving a meaning to 
geodesic distances and geodesic curves in such a setting. 

Yet, the optimal traffic measure Q that one can build from a has the property that iq = \a\ (see 
[1]) and this allows, provided a G , to set the geodesic problem in a usual continuous framework. 
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